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Abstract 

We study whether the neutron skin thickness Ar np of 208 Pb origi- 
nates from the bulk or from the surface of the neutron and proton density 
distributions in mean field models. We find that the size of the bulk 
contribution to Ar np of 208 Pb strongly depends on the slope of the nu- 
clear symmetry energy, while the surface contribution does not. We note 
that most mean field models predict a neutron density for 208 Pb between 
the halo and skin type limits. We investigate the dependence of parity- 
violating electron scattering at the kinematics of the PREX experiment 
on the shape of the nucleon densities predicted by the mean field mod- 
els for 208 Pb. We find an approximate formula for the parity-violating 
asymmetry in terms of the central radius and the surface diffuseness of 
the nucleon densities of 208 Pb in these models. 

1 Introduction 

The neutron skin thickness of nuclei is defined as the difference between the 
root mean square radius of neutrons and protons, 

Ar np = (r 2 )^ 2 - (rX /2 ■ (1) 

The extraction of neutron radii and neutron skins from the experiment is in 
general dependent on the shape of the neutron distribution used in the analysis 
[TJ |31 S] . The data typically do not indicate unambiguosly if the difference 
between the peripheral neutron and proton densities is caused by an extended 
bulk radius of the neutron density, by a modification of the width of the surface, 
or by some combination of both effects. In particular, the neutron skin Ar np of 
208 Pb is nowadays attracting significant interest in both experiment and theory 
since it has a close relation with the density dependence of the nuclear symmetry 
energy and with the equation of state of neutron- rich matter [5], which have 
a large impact in diverse problems of nuclear physics and astrophysics where 
neutron- rich matter is involved [5J El E] . 
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2 Method 



The analysis of bulk and surface contributions to the neutron skin thickness of a 
nucleus requires proper definitions of these quantities based on nuclear density 
distributions [5] . Using the standard definitions of the equivalent sharp radius 
R and surface width parameter b of a nucleon density profile, we have shown in 
Refs. [HI 0] that one can write Ar np = Ar^p lk + Ar™ rf , where 



and 




We recall that R q stands for the radius of a uniform sharp distribution whose 
density equals the bulk value of the actual density and that has the same number 
of particles. Therefore, one has a natural splitting of Ar np in terms of a bulk 
part (0) independent of surface properties, and a part §5$ of surface origin [510]. 



3 Results 

3.1 Model predictions for bulk and surface contributions 

Nonrelativistic (NRMF) and relativistic (RMF) mean field models often differ 
in their predictions for properties of asymmetric nuclear matter. A common 
example is the value predicted for the density slope L of the nuclear symmetry 
energy at saturation density, which may show large discrepancies in the MF 
interactions. The L parameter is defined as 



3po 
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dp 



(4) 



where c sym (p) is the symmetry energy and po the nuclear saturation density. 

In Figure [I] we show the linear correlation of Ar np of 208 Pb with the pa- 
rameter L and demonstrate that it mainly arises from the bulk part of Ar np 
within a large and representative set of mean field models of very different 
nature. Relativistic (RMF) models: G2, NLC, NL-SH, TM1, NL-RAl, NL3, 
NL3*, NL-Z, NL1, DD-ME2 and FSUGold; and nonrelativistic (NRMF) models: 
HFB-8, MSk7, D1S, SGII, DIN, Sk-T6, HFB-17, SLy4, SkM*, SkSM*, SkMP, 
Ska, Sk-Rs and Sk-T4. The original references to the different interactions can 
be found in \§\ [TU] for the Skyrme models, QTj for the Gogny models, and 
OJ [H [H 113 US and Ref. [19] in P2 for the RMF models. From this figure, 
one sees that whereas the bulk contribution to the neutron skin thickness of 
208 Pb changes largely among the different mean-field interactions, the surface 
contribution remains confined to a narrow band of values. The shape of the 
neutron density is more uncertain than the proton density in 208 Pb, and even if 
the neutron rms radius is determined (e.g. in PREX [H]), it can correspond to 
different shapes of the neutron density. As discussed in [4], from the study of 
the two-parameter Fermi functions fitted to the self-consistent densities of the 
MF models, we find that for 208 Pb the central radii C„ vary within the windows 
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Figure 1: Linear correlation of Ar np of 208 Pb with the density slope of the 
nuclear symmetry energy L. The dependence on L of the bulk and surface 
contributions defined in Eqs. © and © is also displayed. 

C„ w 6.7 - 6.85 fm in NRMF and 6.8 - 7 fm in RMF, C p « 6.65 - 6.7 fm in 
NRMF and 6.7 — 6.77 fm in RMF, and that the diffuseness parameters a q vary 
within the windows a n w 0.53 - 0.55 fm in NRMF and 0.55 - 0.59 fm in RMF 
and a p w 0.43 — 0.47 fm. From these results, we can conclude that most of the 
MF models predict a neutron distribution for 208 Pb between the halo-type limit 
(where C n — C p — 0) and the skin-type limit (where a n — a p = 0). The halo- type 
is supported by models with a very soft symmetry energy. Models with a stiffer 
symmetry energy (larger L values) have C n — C p more and more different from 
zero. However, a pure skin-type is not found in any mean-field model as a n — a p 
is always non-vanishing. 

3.2 Parity violating electron scattering at the kinematics 
of the PREX experiment 

Parity-violating electron scattering (PVES) probes the neutron density in a nu- 
cleus since the Z° boson couples mainly to neutrons [T7j . The PREX experiment 
at JLab [T5] aims to determine the neutron rms radius of 208 Pb by PVES. We 
have investigated the dependence of PVES on the shape of the neutron and 
proton densities of 208 Pb within the MF models. To this end, we compute the 
parity- violating asymmetry 

/ da+ der_ \ I ( d<r+ do- \ 

Alr = \ln~ln)/\-dn + ln) (5) 

at the PREX kinematics [T8] . We obtain the differentical cross sections dcr± j d£l 
for right and left-handed electrons by performing the exact phase shift analysis 
(DBWA) of the Dirac equation [H] for incident electrons moving in the po- 
tentials V±(r) = Vbouiomb ( r ) ± V wea k(r) [3]. See Ref. [I] for the details of the 
calculations. We display in Fig. [2] the results for Alr as a function of C n — C p 
and a n — a p of the two parameter Fermi distributions fitted to the 208 Pb MF 
densities. We have found that the results can be reasonably parametrized by 
the formula 10 7 Alr ~ a + f3(C n — C p ) + j(a n — a p ) with a — 7.33, (3 = —2.45 
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Figure 2: Parity violating asymmetry (DWBA) in 208 Pb for 1 GeV electrons 
at 5° scattering angle. MF results and those obtained with the parametrized 
formula given in the text. 

fm, -1 and 7 = —3.62 fin -1 , which is depicted by the crosses in the figure. Re- 
cently, in [50] we have analyzed systematically the correlations of Alr with the 
neutron skin thickness of 208 Pb and with the slope of the symmetry energy in 
the nuclear MF models. 



4 Conclusions 

We have found that the known linear correlation of Ar np of 208 Pb with the den- 
sity derivative of the nuclear symmetry energy arises from the bulk part of Ar np . 
This implies that an experimental determination of R n of 208 Pb could be as use- 
ful for constraining L as a determination of (r 2 )}/ 2 . MF models can accomodate 
the halo-type distribution in 208 Pb if the symmetry energy is very soft but do 
not support a purely skin-type distribution. We find a simple parametrization 
of Alr in terms of C„ — C p and a n — a p that would provide a new correla- 
tion between the central radius and the surface diffuseness of the distribution 
of neutrons in 208 Pb assuming the proton density known from experiment. 
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